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ABSTRACT 

We present a re-analysis of the European Large Area ISO Survey (ELAIS) 90 /im 
observations carried out with ISOPHOT, an instrument on board the ESA's Infrared 
Space Observatory (ISO). With more than 12deg 2 , the ELAIS survey is the largest 
area covered by ISO in a single program and is about one order of magnitude deeper 
than the IRAS 100 fim survey. The data analysis is presented and was mainly per- 
formed with the Phot Interactive Analysis software (PIA, Gabriel et al. 1997) but 
using the pairwise method of Stickel et al. (2003) for signal processing from ERD 
(Edited Raw Data) to SCP (Signal per Chopper Plateau). The ELAIS 90 ^m cata- 
logue contains 229 reliable sources with fluxes larger than 70 mJy and is available at 
http://www.blackwell-synergy.com Number counts are presented and show an excess 
above the no-evolution model prediction. This confirms the strong evolution detected 
at short er(15/xm) and longer (170 /im) wavelengths in other ISO surveys. The ELAIS 
counts are in agreement with previous works at 90 fira and in particular with the 
deeper counts extracted from the Lockman hole observations (Rodighiero et al. 2003). 
Comparison with recent evolutionary models show that the models of Franceschini et 
al. (2001) and Guiderdoni et al. (1998) which includes a heavily-extinguished popu- 
lation of galaxies give the best fit to the data. Deeper observations are nevertheless 
required to better discriminate between the model predictions in the far-infrared and 
are scheduled with the Spitzer Space Telescope (e.g. Lonsdale et al. 2003) which al- 
ready started operating and will also be performed by ASTRO-F (e.g. Pearson et al. 
2004). 
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1 INTRODUCTION 

Strong evolution has been detected in the infrared regime 
based on IRAS number counts at 12, 25, 60 and 100 )im 
(Hacking & Houck 1987, Hacking, Condon & Houck 1987, 
Hacking & Soifer 1991 , Oliver, Rowan-Robinson & Saunders 
1992, Bertin, Dennefeld & Moshir 1997) which show an ex- 
cess of galaxies compared to the no-evolution scenario. These 
findings were recently confirmed with much deeper surveys 
carried out with the ISOPHOT instrument on-board the 
Infrared Space Observatory {ISO) (Kessler et al. 1996) at 
90 and 170 /im (Kawara et al. 1998, Puget et al. 1999, Efs- 
tathiou et al. 2000, Linden- V0rnle et al. 2000, Juvela, Mat- 
tila & Lemke 2000, Matsuhara et al. 2000, Dole et al. 2001). 
ISO also detected a substantial number of faint sources, con- 
sistent with strong evolution from 15 /im number counts (see 
e.g. Elbaz et al. 1999, Gruppioni et al. 2002). Differential 
counts obtained from several independent 15 /im ISOCAM 
surveys show a remarkable upturn at Si 5 < 3 mjy and an ex- 
cess of a factor 10 at the faintest flux above the no-evolution 
predictions. 

In addition to the excess of galaxies detected by ISO 
surveys from the mid-infrared to the far-infrared, the ob- 
servational constraints set by the discovery of the cosmic 
infrared background (CIB) (see Hauser & Dwek 2001 for 
a review and references therein) together with deep sub- 
millimetre surveys (Hughes and Dunlop 1998, Barger et al. 
1998, Eales et al. 2000, Scott et al. 2002, Webb et al. 2003) 
are dramatically increasing the development of new scenar- 
ios of galaxy formation and evolution (Pearson & Rowan- 
Robinson 1996, Guiderdoni et al. 1998, Devriendt & Guider- 
doni 2000, Rowan-Robinson 2001, Franceschini et al. 2001, 
Takeuchi et al. 2001, Pearson 2001, Wang 2002, Lagache et 
al. 2003, Xu et al. 2003). 

The ELAIS survey (for an overview see Oliver et al. 
2000, Paper I) was the largest open time project conducted 
by ISO. This survey consists of more than 12 deg 2 of the 
sky surveyed at 15 and 90 /im, nearly 6 deg 2 at 6.7 /im and 
1 deg 2 at 175 /im (i.e. the FIRBACK survey, see Puget et al. 
1999) in four high Ecliptic latitude (\f3\ >40°) regions with 
low IRAS 100pm sky brightness (< 1.5 MJy sr" 1 ). In this 
work, we present a 90 /im analysis and source counts limited 
to the four large areas, three in the northern hemisphere 
(Nl, N2 and N3), and one in the southern hemisphere (SI). 
Preliminary results of the ELAIS survey at 90 /im based on 
the Quick Look analysis and the brightest sources were pre- 
sented in Efstathiou et al. 2000 (hereafter referred to as Pa- 
per III). 

The paper is organized as follows: in Section [5] we de- 
scribe the observations and the data reduction based on the 
analysis of the distribution of consecutive read-outs of the 
detector instead of using the whole ramp. After the source 
extraction (Section we search in vain for solar system 
objects to remove them from the source list and since some 
could be useful for calibration purposes. In Section^] we esti- 
mate the completeness of the survey, source flux and position 
accuracies and the Eddington bias correction from Monte- 
Carlo simulations of artificial sources on the final maps. The 
final catalog of sources is presented in Section^ We compare 
the ISOPHOT calibration for all standard stars observed at 
90 /im with model predictions and for the sources detected 
in the survey with IRAS values (Section Temperatures 



from colour ratios between 90 and 170 /im for sources also 
detected in the FIRBACK survey are computed in Section 
Q After computing the structure noise (Section |5J in the 
ELAIS fields we present number counts (Section [5J which 
are compared with other works at 90 /im and to evolution- 
ary models before the summary and discussion of our results 
in section [TUl 

2 OBSERVATIONS AND DATA PROCESSING 

2.1 Observations 

The 90 /xm ELAIS data consist of 13 to 20 P22 staring raster 
maps performed with the 3x3 array detector C100 of the 
ISOPHOT instrument (Lemke et al. 1996, for an overview 
see the ISOPHOT Handbook by Laureijs et al. 2003) on 
board ISO. The pixel size on the sky of the C100 detector 
is 43". 5 x 43". 5 and the distance between the pixel centers 
are 46". The ISOPHOT filter-band C 90 with a reference 
wavelength of 90 /im and a width of 51 /tm was used. At 
this wavelength, the full width at half maximum (FWHM) 
of the beam profile is 50 arcsec. Each raster map covers 
typically 20x40 arcmin 2 . Table Al in Paper I provides full 
details on the observations. The Nl, N2, N3 and SI fields 
cover 2.74, 2.98, 2.16 and 4. 15 deg 2 on the sky respectively 
i.e. 12.03 deg 2 in total. The exposure time was 20s but a 
number of sub-fields (representing about 17% of the whole 
survey area) were re-observed with 12s exposures (see Fig. 
14 in Paper I for the survey coverage). 

2.2 Signal processing 

The data were first processed with the Phot Interactive 
Analysis (PIA) software (Gabriel et al. 1997) version 9.1 
using the OLP10 calibration files modified by the inclusion 
of the new dark signal correction (del Burgo et al. 2003a, 
2003b). The data reduction from ERD (Edited Raw Data) 
to SCP (Signal per Chopper Plateau) was performed using 
the pairwise method of Stickel et al. (2003) which was also 
used by Juvela et al. (2000). The signal derived from the dis- 
tribution of the difference between consecutive read-outs is 
used instead of making linear fits to the whole ramps. After 
rejecting the first 10 per cent of the data stream which may 
be affected by transient, the unweighted myriad technique 
(Kalluri & Arce 1998) was used as a robust estimator of 
the pairwise distribution for each raster position. The dis- 
tribution was assumed to be Cauchy (a type of a-stable 
distribution like Gaussians but with a heavier-tailed distri- 
bution) to take into account the presence of glitches in the 
tail of the distribution. 

2.3 Calibration 

Each raster was preceded and followed by an FCS (Faint 
Calibration Source) measurement. However, the calibra- 
tion of the on-sky measurements was made using the sec- 
ond FCS only, performed immediately after the raster and 
with a power chosen to reproduce the intensity of the sky 
background of the measurement. The second FCS gener- 
ally shows a smaller transient behavior than the first one, 
providing a more accurate measurement. For each field, the 
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relative uncertainty coming from the FCS calibration was 
computed as the mean absolute deviation of the average sky 
background of all rasters, which results to be 7 per cent. 



2.4 Flat-fielding and Mapping 

Differences of up to 20% in the overall levels of the data 
streams of the detector pixels were noticed after the flux cal- 
ibration. This behaviour is most likely resulting from pixel- 
to- pixel sensitivity differences, which moreover appeared to 
be time-dependent. 

To correct for this, the pixel data streams were slightly 
smoothed and filtered to remove sources. At each raster 
point the mean of the filtered pixel values was computed. 
The sequence of ratios of the mean and the individual pixel 
value at each raster point was fitted with a robust polyno- 
mial to give the smooth correction function for each pixel. If 
remaining time trends were still noticeable after correcting 
the individual pixels to the common mean (by multiplica- 
tion), the procedure was repeated but the filtered data val- 
ues from all detector pixels were simultaneously fitted with 
a robust low order polynomial. This removes any time trend 
still present after rescaling the pixel data streams to the 
common mean. 

This combined method is highly effective in removing 
pixel-to-pixel sensitivity differences and time trends in the 
data streams. 

The whole field map was built from the Jy /pixel values 
for each field using the drizzle mapping method under IRAF 
(Tody 1993) with a pixel size of 30 arcsec and the default 
shrink factor (0.65). The drizzle method (Fruchter & Hook 
1997, 2002) allows to consider the exact size of pixels and 
gaps between them (see Sect. 12. ll . 



3 SOURCE DETECTION 

The source detection was performed using the SExtractor 
software version 2.2.2 (Bertin & Arnouts 1996) on the final 
maps. The sky background was computed in a grid of 15 x 15 
square pixel (i.e. 7.5 x 7.5 square arc-minute) and SExtractor 
was run with a detection threshold of 1.8 a and a minimum 
number of pixel equal to 2. The flux in a circular aperture 
(FLUX_APER) of 6 pixels (i.e. 180 arcsec) diameter was 
used. 



3.1 Search for solar system objects in the ELAIS 
fields 

Although the ELAIS fields are at high ecliptic latitude (for 
low zodiacal background), there could still be objects from 
the solar system inside the fields. The search for these tar- 
gets had two purposes: 1. Cleaning of the ELAIS source 
list from moving solar system targets. 2. Finding additional 
targets which might be later on used for independent flux 
calibration purposes (Muller et al. 2002; Miiller & Lagerros 
1998, 2002). As the raster maps were observed at different 
periods during the ISO mission, we used the exact date and 
time at which they were obtained to search the databases of 



the Minor Planet Center 1 . For our search we included more 
than 150 000 asteroids with reliable orbital elements (num- 
bered asteroids and unnumbered, multi-apparition objects) , 
more than 200 comets and the planets and their satellites. 
A search radius which was slightly larger than the actual 
ELAIS fields was used to account for the geocentric to ISO- 
centric parallax errors (the position calculations were done 
in the geocentric frame) . A geocentric to /SOcentric parallax 
of 10 arcmin covers all asteroids beyond 0.15 AU from Earth, 
i.e. more than 99% of all known asteroids, but we allowed 
for parallaxes of up to 30 arcmin to also account for possible 
ephemeris uncertainties and asteroid movements during the 
observations. This means that the 2 deg search radius for 
each ELAIS field included a very large safety margin. 

One asteroid and one comet were selected to be possi- 
bly seen in one of the ELAIS rasters. We computed their 
movements during the observation: 12 arcsec and 1.5 ar- 
cmin. The ISO parallax correction was less than 1 arcmin in 
both cases. The two objects were finally found to be outside 
the ELAIS field when we repeated the ephemeris calculation 
with a more sophisticated N-body tool in the /SOcentric 
frame. From our analysis we concluded that there are no 
known solar system sources in the ELAIS 90 /im data. 



4 COMPLETENESS, SOURCE FLUX 
UNCERTAINTY AND POSITION 
ACCURACY 

To estimate several quantities such as completeness, flux 
and positional uncertainties we adopt a similar approach as 
Dole et al. (2001) for FIRBACK based on the addition of 
artificial sources to the data. Artificial sources were added 
at random positions on the final map of each field using 
the 90 /im theoretical footprint scaled by a certain factor to 
simulate sources with a known flux. 

In practice, to keep only the noise on the im- 
ages, objects detected by SExtractor were first removed 
from the images for the simulations (subtracting the im- 
age obtained with the SExtractor option CHECKIM- 
AGE_TYPE=OB JECTS) . The source position can fall any- 
where on a pixel and the pixelised footprint was computed 
in a square of 5 x 5 pixels providing a spatial extension of 2.5 
x 2.5 arcmin 2 for each source and representing 96 per cent 
of the total flux contained in the theoretical footprint. Sim- 
ulations of 15 artificial sources and the extraction with SEx- 
tractor using the same parameters as for the survey sources 
were repeated 300 times for each field giving a total of 4500 
simulated sources at each flux level equal to 50, 60, 70, 80, 
90, 100, 125, 150, 200, 300, 400, 500, 750 and 1000 mjy. 

4.1 Positional accuracy 

The positional accuracy can be estimated from the statisti- 
cal analysis of the distances between recovered sources and 
the exact position of simulated sources. Figure shows the 
histogrammes of distances between extracted and simulated 
sources with fluxes equal to 100, 200 mjy and 500 mjy for 
the N2 field. The peak of the distribution is around 8 arcsec 

1 http:/ /cfa-www. harvard.edu/cfa/ps/mpc. html 
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Figure 1. Distribution of the distance between recovered and 
simulated sources with a flux equal to 100, 200 and 500 mjy 
in plots 'a','b' and 'c' respectively in N2. The mean distance is 
indicated as a dashed line in each graph and is equal to 8.4, 5.4 
and 3.3 arcsec for 100, 200 and 500 mjy, respectively. 



Figure 2. Histogrammes of the recovered to input flux for simu- 
lated sources with a flux equal to 100, 200 and 500 mjy in plots 
'a','b' and 'c' respectively in N2. The mean unity is indicated as 
a dashed line. Solid lines are Gaussian fits to the distribution. 



for 100 mjy sources and below 5 arcsec for sources brighter 
than 200 mjy. 

The absolute pointing error of ISO represents only a 
small additional uncertainty as it was better than a few arc- 
sec all along the mission (Kessler 2000). 

4.2 Flux uncertainties 

Histogrammes of the recovered to input flux of simulated 
sources are shown on Fig.|2|for N2 at 100, 200 and 500 mjy. 
At each flux level, the distribution was fitted with a Gaus- 
sian whose a gives an estimate of the photometric accuracy. 
Fig. gives the variation of a as function of flux level de- 



rived from simulated sources detected with SExtractor with 
a signal-to-noise ^ 3. The uncertainty on the recovered flux 
is typically 30 per cent at 100 mjy and decreases to less 
than 10 per cent above 400 mjy. 



4.3 Completeness 

The completeness of the survey is computed as the ratio of 
the number of recovered sources with signal-to-noise ratio 
above 3 to the total number of simulated sources and is 
shown on Fig. Eh for N2. The completeness is almost 100 
per cent down to 150 mjy and decreases to 77 per cent at 
100 mjy. 
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Figure 3. a: Completeness of the survey as function of flux (in 
logarithmic scale) derived from simulation of artificial sources in 
N2. The Completeness was computed for sources with signal-to- 
noise ratio better than 3. ft: Flux uncertainty derived from Gaus- 
sian fitting of the distribution of recovered to input flux as func- 
tion of input flux. 
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Figure 4. The Eddington bias as function of flux (in logarithmic 
scale) for the ELAIS survey computed as the centre of the gaus- 
sian fit to the distribution of measured to input source flux (see 
Sect. 14751 . The mean correction (diamonds) is the average of the 
values for the four ELAIS fields (plus signs). The solid line is the 
result of a polynomial fit of degree 4. Error bars are the standard 
deviations at each flux level. 



mjy, respectively. The correction for the Eddington bias was 
directly performed on the source flux (while the usual way 
is to correct the number counts assuming a certain power 
law (see references above)). This provides corrected source 
catalogues and does not need any assumption on the distri- 
bution of source flux to apply the correction to the number 
counts. 



4.4 Eddington bias 

Noise on the images is responsible for an excess in the num- 
ber counts as it will create an overestimate of fluxes. This 
effect, know as Eddington bias (Eddington 1913), is similar 
to Malmquist bias, which refers to fluctuation in intrinsic 
rather than measured quantities (see e.g. Teerikorpi 1998). 

The proper determination of the bias plays an impor- 
tant role in the estimation of source flux, the computation 
of number counts and therefore the determination of the 
strength of the evolution seen in the counts as the correction 
dramatically increases towards the faint end of the sample. 

One can estimate the Eddington bias analytically as- 
suming a certain power-law and adding an appropriate flux 
dispersion like in e.g. Murdoch et al. (1973) for an under- 
lying Euclidean slope, Oliver et al. (1995) and Dole et al. 
(2001) who all assumed Gaussian noise. One can also use 
a Monte-Carlo approach like Bertin, Dennefeld & Moshir 
(1997). 

A more realistic estimate of the bias can be obtained 
from simulations performed on the maps themselves to es- 
timate the correction. A mean correction of the Eddington 
bias was computed for the four fields and is presented in 
Fig. 0] as a polynomial fit to the centres of the Gaussian fits 
to the distributions of measured to input flux (Sect 14.21 and 
Fig.0. The bias is less than 34 and 13% above 70 and 100 



5 ELAIS 90 MICRON FINAL CATALOGUE 

To check the reliability of the sources detected with SEx- 
tractor, the classification presented in Paper III was used to 
check the reliability of the detected sources. Five persons 
eyeballed all the detections above 1.5 a of the sky back- 
ground detected along the data streams of individual pix- 
els. Only sources with detections classified at least twice 
as probable sources within a circle of 150 arcsec radius, a 
signal-to-noise ^ 3 and a flux 70 mjy (i.e. above the 3-cr 
noise level computed in Sect. [HJ were retained for the final 
source list. The selection based on the eyeball classification 
ensures that there are no or few fake sources in our sample. 

The final ELAIS 90 /im source list contains 229 sources 
while the 163 most reliable sources detected in the prelim- 
inary analysis were presented in Paper I with flux uncer- 
tainties estimated to be 40%. The comparison of source 
flux from the final and preliminary analysis is presented 
in Rowan- Robinson et al. (2004) for ISOPHOT and ISO- 
CAM and shows a good agreement. Table Q gives right 
ascension, declination, flux and flux uncertainty (which 
contains the uncertainty given by SExtractor and the er- 
ror coming from the Eddington bias correction) for each 
source. The full version is available in electronic format at 
|http:/ /ww w . blackwell-synerg y.com| 
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Table 1. ELAIS 90 fim source list. Columns give the name of the source using the ELAISP90.JHHMMSS+DDMMSS format according 
to the acronym in the IAU Registry, the right ascension and declination, the flux and flux uncertainty in mjy. The full version of this 
table is available in electronic format at http://www.blackwell-synergy.com 
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00 


33 


18.1 


-44 


08 


28.9 


175 


29 


ELAISP90. 


.J003321-432700 


00 


33 


21.9 


-43 


27 


0.3 


260 


17 


ELAISP90. 


.J003349-441903 


00 


33 


49.3 


-44 


19 


3.7 


79 


18 


ELAISP90. 


J003359-441108 


00 


33 


59.5 


-44 


11 


8.3 


177 


24 


ELAISP90. 


J003415-423205 


00 


34 


15.2 


-42 


32 


5.1 


77 


25 



6 CALIBRATION COMPARISONS 

To check the quality of the calibration at the low sur- 
face brightness level of the ELAIS fields, we compare the 
ISOPHOT calibration with theoretical predictions for stan- 
dard stars (Sect. IoTTi . and with IRAS (Sect. HHfll flux esti- 
mates. 



6.1 Standard stars 

In order to better determine the ELAIS calibration (as well 
as the general ISOPHOT calibration) three stars (HR6132, 
HR6464 and HR5981) close to the ELAIS fields were ob- 
served in mini-raster mode (a 3 x 3 raster with the star 
positioned at the centre of a different pixel in each point- 
ing). The faintest of the stars (HR5981) was observed twice 
on the same ISO orbit. 

To increase the sample of measurements and thus the 
reliability of the comparison, all other standard stars ob- 
served in mini-raster mode at 90 pm were retrieved from the 
ISOPHOT archive. The comparison with two model predic- 
tions was performed. Hammersley et al. (1998) models were 
constructed by fitting near-IR observations performed with 
the Infrared Telescope Facility (IRTF). Cohen et al. (1999) 
constructed empirical stellar spectra in the near and mid- 
infrared based on observations taken from the ground, the 
Kuiper Airborn Observatory, and the IRAS Low Resolution 



Spectrometer. Both predictions were extrapolated to longer 
wavelengths as v 2 . 

Table |5|shows the list of stars and the characteristics of 
the ISOPHOT measurements as well as model predictions. 
The predicted stellar fluxes lie in the range between 60 mjy 
and lOJy at 90 /im. Uncertainties on the models estimates 
are typically 3 and 5% for Hammersley et al. (1998) and 
Cohen et al. (1999), respectively. 

The integration time per pointing in these mini-rasters 
(from 40 to 72 s) is longer than that used for the bulk of the 
ELAIS survey in order to obtain an accurate determination 
of fluxes to establish the ISOPHOT calibration. 

The observations of calibration stars were processed in 
the same way as the survey rasters. The application of a 
method based on celestial standards ultimately depends on 
the accuracy with which the background can be estimated 
and on the accuracy of the fluxes of the sources used as 
calibrators. However, in the case of the small rasters maps 
performed on standard stars, SExtractor fails to compute a 
reliable sky (and therefore star flux) estimate. 

The best way to extract both star and sky background 
estimates for point sources observed in mini-raster mode 
was found to be to use the individual pixel values at each 
raster position weighted by the point-spread function frac- 
tions derived for the ISOPHOT C100 (Laureijs 1999, Moor, 
in preparation). The fraction of the point-spread function 
falling on a C100 pixel situated at a distance d from the 
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1.5- 



point source centre (/ ps f(d)) has been determined for each 
filter at a number of typical distances. The method assumes 2.5 
the source is point-like (the / ps f factors have to be modified 
if the source is extended) and centered on the detector pixels 
which is the case for standard stars. ISOPHOT values were 2.0 - 

color-corrected according to the spectral type of these stars. 

Results of the comparison are given in Table |5] as the 
ratio between measured (based on the FCS) and theoretical 
fluxes. When two model predictions were available, we used 
their weighted mean to compute the ratio. The measurement 
of a star (HR7451) with a very low predicted flux (7.5 mjy) 
was excluded from the comparison. The two measurements 
of the brightest star (HR5340) are in very good agreement. 

ISOPHOT fluxes are on average higher than the pre- 
dicted ones and the weighted mean ratio is 1.06±0.02. It is 
unclear whether this discrepancy is coming from the differ- 
ences in the observing setup used for standard stars and the 
ELAIS survey. The difference between the FCS calibration 
and the model prediction for stars is shown on Figs. 8, 9 and qq q g ^ q -| 5 2 25 

10. IRAS Flux (Jy) 




6.2 Comparison with IRAS sources 

While the ELAIS fields were chosen to avoid strong 12 )im 
infrared sources, there are a number of IRAS 100 fim sources 
detected in the survey. All common sources have low (the 
flux is an upper limit) or intermediate IRAS quality flags 
(Moshir, Kopman, & Conrow 1992). Figure [£] shows the 
comparison with the Faint Sourve Catalogue (FSC) which 
is known to be more accurate than the Point Source Cata- 
logue at faint level. Only sources with quality flags equal to 
2 (intermediate accuracy) were selected for the comparison 
and this represents 21 IRAS sources. 

Colour correction factors were computed from the 
IRAS 4-band Spectral Energy Distribution and IRAS and 
ISOPHOT filter profiles. 

The mean ratio of ISOPHOT to IRAS Hux is 0.76 with 
a standard deviation of 0.17 and shows a discrepancy with 
the model prediction comparison (Sect. l6~Tt . 

However it should be noted that there is a tendency for 
the IRAS FSC to overestimate fluxes near the FSC thresh- 
old at 60 and 100 ^m (Moshir et al 1992). Since all the 
IRA S sources in Fig. E| have 5(100) < 3 Jy (and most have 
5(100) < Uy), this would be sufficient to explain the dis- 
crepancy noted above. 



7 CORRELATION WITH FIRBACK 

We looked for FIRBACK identifications (Dole et al. 2001) 
of our 90 jim source sample within a circle of 188 arcsec ra- 
dius (i.e. sqrt(2) x (89.4+43.5) where 89.4 and 43.5 are the 
pixel size of the C100 and C200 detectors, respectively) , 
also including the complementary FIRBACK source cata- 
log which provides sources with fluxes down to 135 mjy. If 
several sources were selected, the closest identification was 
used. In the common area of the two surveys, 53 out of 102 
and 21 out of 55 FIRBACK sources were identified at 90 fim 
in the Nl and N2 fields, respectively. 

Since ISOPHOT fluxes refer to a spectrum with vF v — 
constant, colour temperature Tc were computed correcting 



Figure 5. Comparison of ISOPHOT and IRASFSC fluxes at 
90 micron for 21 common sources with intermediate IRAS qual- 
ity flags (i.e. 2). The 90 fim fluxes of the IRAS sources are esti- 
mated by linearly interpolating in linear space between the colour- 
corrected 60 and 100 fim fluxes. 



20 



15 



5 



o i , . , m 

10 15 20 25 30 

170/90 micron Colour Temperature [K] 

Figure 6. The 170fim/90fim colour temperature distribution 
for 53 and 21 sources detected at both wavelengths in Nl and N2 
fields respectively. An emissivity index (3 = 2 was adopted. 



the 90 and 170 fim fluxes in the two band-passes for a mod- 
ified blackbody function with an emissivity index j3 = 2. 

The resulting distribution of colour temperatures 
(Fig. |SJ is centered around Tc = 19 K with most of sources 
lying in the range 15-25K. These values are in favour of 
the presence of a cold component in low redshift galaxies 
(Rowan-Robinson et al. (2004) derived a median redshift 
of 0.15 for the 90 pm sample) which was first detected in 
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Table 2. The list of standard stars used to check the FCS calibration with theoretical values. Colums are : the target dedicated time 
(TDT) number of measurement, name of stars, exposure time, size of the mini-rasters in steps of 46 arcsec, model predictions from Cohen 
et al. 1998 (-Fmc) an d Hammersley et al. 1998 (-Fpjr) and ISOPHOT measurements (-Fphot) ar >d their respective uncertainties (eMC i 
epjj and ep hot ) are indicated. "Ratio" is the ratio of measured to predicted fluxes. When two model predictions were available, we used 
their weighted mean to compute the ratio. Errors on the ratio are given in eR at ; . The weighted mean ratio is 1.06±0.02. 













Models 




Measurements 






Measurement 


Name 


Exposure 


Size 




e MC 


^PH 


epH 


-Fphot 


ephot 


Ratio 


e Ratio 


TDT number 




sec 




Jy 


Jy 


Jy 


Jy 


Jy 


Jy 




% 


08602417 


HR5340 


37.00 


5 


9.303 


0.528 


9.029 


0.300 


9.54 


0.34 


1.05 


0.05 


10503417 


HR6705 


72.00 


5 


2.012 


0.115 


1.904 


0.066 


2.02 


0.14 


1.05 


0.08 


27502117 


HR5340 


72.00 


5 


9.303 


0.528 


9.029 


0.300 


9.52 


0.23 


1.05 


0.04 


29301005 


HR7310 


72.00 


5 


0.258 


0.015 


0.268 


0.009 


0.33 


0.02 


1.24 


0.08 


39103002 


HR8775 


72.00 


5 


4.957 


0.282 


5.096 


0.184 


5.49 


0.15 


1.09 


0.04 


65701318 


HR1654 


72.00 


3 


0.713 


0.042 






0.74 


0.02 


1.04 


0.07 


72701418 


HR7980 


72.00 


3 


0.517 


0.031 






0.48 


0.02 


0.93 


0.07 


77200361 


HR5981 


40.00 


3 






0.063 


0.002 


0.07 


0.02 


1.11 


0.32 


77200364 


HR5981 


40.00 


3 






0.063 


0.002 


0.07 


0.01 


1.11 


0.16 


78300465 


HR6464 


40.00 


3 






0.120 


0.004 


0.13 


0.02 


1.08 


0.17 


78300677 


HR6132 


40.00 


3 






0.288 


0.001 


0.32 


0.03 


1.11 


0.10 



the ISO Serendipity Survey (Stickel et al. 1998, 2000, 2001). 
Dunne and Eales (2001) also recently measured a cold com- 
ponent (20-21K) in their sample of 17 galaxies combining 
IRAS and SCUBA observations at 450 /im. This is consistent 
with the analysis of COBE / DIRBE data for the MW Galaxy 
(Sodroski et al. 1994). A more detailed analysis with spectral 
energy distribution fitting of the ELAIS sources from the op- 
tical to the FIR domain is presented in Rowan- Robinson et 
al. (2004). 



8 SOURCE CONFUSION ESTIMATES 

Estimates of the confusion noise relies on the direct measure- 
ment of structure noise, Af s t r . The structure noise is calcu- 
lated via the so-called structure function S, which measures 
the average brightness fluctuations for a specific measure- 
ment configuration see e.g. Herbstmeier et al. (1998) : 



S* k (0) = 



(1) 



where B(x) is the measured sky brightness at the position 
'x', 9 is the separation between the target and reference 
positions, 'k' is the number of reference positions, # r s are 
the vectors to the reference positions relative to the target. 
The average is taken over the whole map. The actual val- 
ues of ; -s are determined by the geometry of the measure- 
ment configuration. The structure noise is calculated from 
the structure function: 



A str = ^/s^ x n 



(2) 



where Q is the effective solid angle of the aperture. The 
structure noise contains the contribution of the sky bright- 
ness fluctuations (confusion noise, iV conl ) and that of the 
average instrument noise Ainst. As was shown by Kiss et 
al. (2001) the relation between these quantities can be well 
approximated by the following formula for ISOPHOT mea- 
surements: 




NLa + 2Afn St 



(3) 



40 60 
N str (mJy) 

Figure 7. Distribution of the structure noise in the ELAIS Nl 
field. The solid line (histogram) represents the measured Nstr val- 
ues over the whole field. The dashed line represents the Gaussians 
fitted to the low amplitude regime (see text for details). 



We derived the distribution of structure noise for individual 
pixel pairs (without averaging in space in Eq. for the 
ELAIS fields Nl, N2, N3 and SI. The results of the ELAIS 
Nl field are presented on Fig. (the N2, N3 and SI fields 
have a similar distribution). 

As seen in this figure the distribution of iV s t r is 
Gaussian- like, with an extended tail toward high structure 
noise values. In brighter cirrus regions the distribution of cir- 
rus fluctuations can be well separated from the fluctuation of 
the CFIRB and the contribution of the instrument noise. In 
those fields the fluctuation distribution can be well described 
by a Gaussian one (Kiss et al. 2003). In the following we as- 
sume that this distribution can be applied for the fluctutions 
in the faint ELAIS fields as well. Fig. [7] can be used to esti- 
mate the source confusion limits of the ELAIS fields. We fit- 
ted a Gaussian to the lower noise regime (iV s t r < 20mJy) of 
the iV a tr distribution of the Nl, N2, N3 and SI fields, which 
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Table 3. Integral number counts in the ELAIS survey. S is the 
flux in mjy; N is the number of sources per deg 2 with flux larger 
than S; The upper and lower uncertainties in the counts which 
come from Poisson error and the correction of the Eddington bias 
are indicated. N_s is the number of sources with flux larger than 
S. 



s 


N(> 5) 


N_s 


mjy 


deg - 2 




95 


17.02±4.50/4.30 


185 


177 


7.93±1.19/1.83 


90 


330 


2.70±0.74/0.74 


30 


614 


0.77±0.39/0.39 


9 



Table 4. Normalized differential number counts. Columns are : 
flux range for each bin in mjy; the bin centre; dN/dS X S 2 5 
are the number counts per bin normalized to the Euclidian law in 
deg - 2 Jy 1 ' 5 (the bin centre in linear scale was used for the normal- 
isation); The upper and lower uncertainties in the counts which 
come from Poisson error and the correction of the Eddington bias 
are indicated. N_s is the number of sources per bin. 



flux bin 


bin centre 


dN/dS x S 2 ' 5 


Nj3 


mjy 


mjy 


deg -2 Jy 1,5 




95-176 


135 


0.74±0.20/0.22 


95 


176-329 


253 


l.ll±0.24/0.24 


60 


329-613 


171 


1.04±0.27/0.27 


21 


613-1142 


877 


0.85±0.32/0.32 


7 



resulted in a a of 14.8, 12.8, 13.4, and 17.1 mjy, respec- 
tively. With the point spread function fraction coefficient of 
/psf = 0.61 for a C100 camera pixel at 90 /im the 3— a source 
confusion limit is ~70mJy. Moreover, it is worth to mention 
that the cosmic far-infrared background has an expected 
fluctuation power of ~7mJy at this wavelength for the C100 
camera detector pixels (Kiss et al. 2001), which contributes 
to the final width of the N str distribution. Eliminating this 
value from the width of the Gaussians, the remaining contri- 
butions of the cirrus fluctuations and the instrument noise 
would be ~60 mjy at 3-<r. 



9 SOURCE COUNTS 
9.1 ELAIS counts 

Integral and normalized differential source counts are given 
in Tables[!|]and[l]for 185 sources brighter than 95mJy (above 
this flux, the uncompleteness and Eddington bias corrections 
are lower than 25 and 15% respectively). Uncertainties in the 
counts represent the contribution of Poisson errors and the 
Eddington bias correction. The possibility that some sources 
could be solar system bodies was rejected in Sec|3] Moreover 
as stated in Paper III, given that there are no bright 12-/im 
sources in the ELAIS fields (to avoid saturating ISOCAM) 
we do not expect any 'photospheric' stars to be detected at 
90 fim as these would have a 90 /im flux ^ 10 mjy. Finally, 
the 3-cr limit of the cirrus and instrumental noise was esti- 
mated in Sect. |S] to be ~60 mjy. Therefore it is very likely 
that all the selected sources above 95 mjy are extragalactic. 
Figure |H| shows ELAIS integral counts (filled circles) 



O) 
CD 
T3 

CO 

A 



° -1 



-1.5 -1.0 -0.5 0.0 0.5 
log[S(90)/Jy] 

Figure 8. Integral source counts at 90 fim for ELAIS (filled cir- 
cles). Error bars are Poisson error plus the uncertainties on the 
Eddington bias. Flux uncertainties are 7% resulting from the use 
of the FCS for calibration purposes (Sect. 2.3). Correction fac- 
tors of 1/1.06 and 1/0.76 derived from the comparison of the 
FCS calibration with standard stars model predictions and with 
IRAS respectively are shown as -0.03 dex and +0.12 dex error 
bars at (-1.3,-2.0). Source counts from Juvela et al. (2000; aster- 
isks), the preliminary analysis of ELAIS (Paper I; open circles), 
Linden- V0rnle et al. (2000; squares), Matsuhara et al. (2000; tri- 
angle) and Rodighiero et al. (2003; diamonds) are shown for com- 
parison. IRAS counts (x symbols) are shown for galaxies in the 
PSCz catalogue. The dashed line is the no-evolution model from 
Franceschini et al. (2001). 



at 90 /im compared with results of Juvela et al. (2000; as- 
terisks), the preliminary analysis of ELAIS (Paper I; open 
circles), Linden- V0rnle et al. (2000; squares), Matsuhara et 
al. (2000) from the Lockman Hole observations (triangle) 
and the new analysis of the Lockman Hole performed by 
Rodighiero et al. (2003) (diamonds). IRAS points (x sym- 
bols) are also shown for galaxies in the PSCz catalogue 
(Saunders et al. 2000) with a selection of galactic latitude 
(|6| > 20) and low IRAS flags (/ qua i < 3) at 100 ^m and 
fluxes brighter than 2Jy as it becomes incomplete at fainter 
level (see Paper III for details). The dashed line is the no- 
evolution model from Franceschini et al. (2001). Correction 
factors of 1/1.06 and 1/0.76 derived from the comparison of 
the FCS calibration with standard stars model predictions 
(Sect. IdTH and with IRAS (Sect. l6~2l are also shown. 

Our results extend IRAS counts by more than one order 
of magnitude. They are in very good agreement with the 
preliminary analysis of ELAIS and confirm the departure 
from the Euclidian slope found in Paper I. Integral counts 
in the range 0.095-lJy are well fitted with a straight line of 
the form: 

log 10 (7V) = (-1.68 ± 0.09) x log 10 (S[Jy]) - (0.43 ± 0.07)(4) 



9.2 Comparison with evolutionary models 

Recent observations in the FIR and sub-millimeter regimes 
have considerably improved evolutionary models in the past 
5 years. In the following, we briefly describe the main charac- 
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Figure 9. Normalized differential source counts for ELAIS (filled 
circles) at 90 (im. Error bars include Poisson error and the uncer- 
tainties on the Eddington bias. Correction factors of 1/1.06 and 
1/0.76 are shown as -0.03 dex and +0.12 dex error bars at (1.5,- 
1.0). Lockman Hole counts from Rodighiero et al. (2003) are also 
shown (diamonds). IRAS counts (crosses) are shown for galaxies 
in the PSCz catalogue. The solid line is the no-evolution model 
from Franceschini et al. (2001). The dashed and dotted line repre- 
sent models A and E of Guiderdoni et al. (1998), respectively. The 
dash-dotted line are the counts predicted by Rowan-Robinson 
(2001) for his cosmological model with Qg = 0.3 and A = 0.7. 



1.5 




-1.5 [ i i i 1 

-2-10 1 2 
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Figure 10. Same as|2]but the dashed-dotted, dashed and dotted 
lines are the models of Pearson (2001), Franceschini et al. (2001) 
and Lagache, Dole & Puget (2003), respectively. 



teristics of evolutionary models by Guiderdoni et al. (1998), 
Rowan- Robinson (2001), Pearson (2001), Franceschini et al. 
(2002), Lagache et al. (2003) and compare the predictions 
to differential number counts measured in the ELAIS survey 
on Figs. El and Ell 

(i) Guiderdoni et al. 1998 have designed a family of semi- 
analytic evolutionary scenarios within the context of hierar- 
chical growth of structures according to the cold dark matter 



model, with prescriptions for dissipative and non-dissipative 
collapses, star formation and feedback. Differences between 
these scenarios only concern the efficiency of star formation 
on a dynamical time-scale, the IMF and the extinction. In 
Figure |5| we compare our results with two of their models: 

Scenario A contains a mix of two broad types of popula- 
tions, one with a 'quiescent' star formation rate, the other 
proceeding in bursts with a high evolution rate and fitting 
the SFR density at low z; 

Scenario E includes an additional population of heavily- 
extinguished galaxies (ULIGs) and is qualified as the best 
fit by Guiderdoni et al. as it nicely reproduces the Cosmic 
Optical Background and the Cosmic Infrared Background. 

If model A is systematically below the measured counts, 
the addition of a ULIGs population shifts the predictions 
upwards and the model E falls in excellent agreement with 
the observations as it was suggested in Paper III based on 
the brightest sources. 

(ii) The models of Rowan-Robinson (2001) include four 
spectral components: infrared cirrus, an M82-like starburst, 
an Arp 220-like starburst, and an active galactic nucleus 
dust torus. The proportion of each spectral type are chosen 
for consistency with IRAS and SCUBA colour-luminosity 
relations and with the fraction of AGNs as a function of 
luminosity in 12 /im samples. 

The prediction of the Rowan- Robinson model for the cos- 
mological model with Qq = 0.3 and A = 0.7 is compared 
with the observed counts on Figure H The model is in good 
agreement with the observations down to fluxes of ~ 200 
mjy within the error bars. At fainter fluxes it gives a slightly 
too high number of sources. 

(iii) The model of Pearson and Rowan- Robinson (1996) 
consists of non-evolving spiral and elliptical components 
mixed with an evolving population of starburst galaxies, ac- 
tive galactic nuclei and a hyperluminous galaxy component. 
The model is in agreement with the source counts at 60 fim 
and the faint radio counts at 1.4 Ghz and provides a good 
estimate of the cosmic infrared background observed with 
COBE at 500 /jm. 

More recently, Pearson (2001) used the framework of the 
Pearson and Rowan-Robinson galaxy evolution model and 
constrains the evolution in the galaxy population with the 
observed counts and background measurement derived from 
ISO and SCUBA observations. Pearson found that a strong 
evolution in both density and luminosity of the ULIG pop- 
ulation can account for the source counts from 15 /im to the 
submillimetre region, as well as explain the peak in cosmic 
infrared background at ~ 140 /im. 

The prediction for this model is also shown on Figure [TU1 
The model provides a good fit to the ELAIS observations 
although it seems to become too high at fluxes fainter 
than ~ 100 mjy compared to the Lockman Hole counts of 
Rodighiero et al. 2003. 

(iv) The model of Franceschini et al. (2001) assumes that 
the extragalactic population is composed of three compo- 
nents with different evolution properties: (1) a non-evolving 
population of spirals; (2) a population of strongly evolving 
starburst galaxies and type-II AGNs; (3) a population of 
type-I AGNs which does not contribute significantly to the 
counts. This model was optimized to reproduce the mid- 
IR counts and redshift distribution. In particular, the two 
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components of the fast evolving population were required to 
reproduce the shape of the 15/im counts. 

The model of Franceschini is plotted on Figure 1101 and 
gives a good estimate of the observed counts. The increase 
in number counts seen at fluxes fainter than ~ 100 mjy is 
probably the far- infrared counterpart of the upturn detected 
in the mid-infrared (Elbaz et al. 1999, Chary and Elbaz 2001, 
Mazzei et al. 2001; Serjeant et al. 2000). The ELAIS 90 (im 
data do not allow us to test if this predicted feature is real or 
not but Rodighiero et al. (2003) have shown it is compatible 
within the error bars with the faint counts in the "Lockman 
Hole". 

(v) Lagache, Dole & Puget (2003) have developed a phe- 
nomenological model which fits all the existing counts and 
redshift distributions from the mid-infrared to the submil- 
limetre range together with the intensity and fluctuation of 
the cosmic infrared background. Their model is based on the 
evolution of galaxy luminosity function with redshift for a 
population of starburts and normal galaxies. 

The model of Lagache, Dole & Puget 2003 (shown on Fig- 
ure EI] is compatible with ELAIS counts (although slightly 
higher) around 1 Jy and with a larger discrepancy at fainter 
level (S ^ 200 mjy) where the model continues to increase 
while the observed counts decrease. 



10 SUMMARY AND DISCUSSION 

We have used a new method to reduce ISOPHOT measure- 
ments in the 4 main areas of the ELAIS survey at 90 /im. 
With a total area of more than 12deg 2 , the ELAIS survey 
represents the largest area covered in a single programme 
with ISO. 

The relative uncertainty in flux coming from the FCS 
calibration estimated from the sky background level differ- 
ences of all rasters is 7 per cent. 

On the one hand, the comparison of measured fluxes 
with models for standard stars shows a strong correlation 
with a mean ratio of ISOPHOT to model values of 1.06 ± 
0.02. On the other hand, the comparison with the IRASFSC 
catalogue for IRAS sources detected in the survey gives a 
mean ratio of ISOPHOT to IRAS values equal to 0.76±0.17. 

Simulations of artificial sources on the final maps span- 
ning a wide range of flux were used to estimate flux and po- 
sitional uncertainties, completeness and the Eddington bias 
corrections. The completeness of the survey is about 80% at 
100 mjy. 

We present a source list of 237 reliable sources with 
fluxes larger than 70 mjy, signal-to-noise 3 for the 4 large 
ELAIS fields. The full version of the catalogue is available 
at http:/ /www. black well- synergy, com 

Sources detected at 90 and 170 /im in the FIRBACK 
survey (Dole et al. 2001) have an average colour temperature 
of Tc = 19 K with all sources lying in the range 13-25K in 
agreement with Stickel et al. (2001) in the ISO Serendipity 
Survey (Stickel et al. 1998). 

The ELAIS counts extend the IRAS counts by more 
than one order of magnitude in flux and show significant 
departure from the no-evolution model as detected in other 
ISO surveys from the mid- to the far-infrared. 

There is in general a good agreement between ELAIS 
and other 90 /im source counts and in particular with the 



deeper counts measured in the Lockman Hole (Rodighiero et 
al. 2003). Differential number counts measured in the ELAIS 
regions at 90 /im are compared to recent evolutionary mod- 
els. Among few models which were compared to our counts, 
the model of Franceschini et al. (2001) and the scenario E 
of Guiderdoni et al. (1998) give the best agreement with the 
observations. 

However, the latter model is a factor of ~ 2.5 below 
the counts measured at 170 /im in two of the ELAIS re- 
gions (Dole et al. 2001) related to the present paper. On the 
other hand Matsuhara et al. (2000) found that the scenario 
E model prediction of Guiderdoni et al. is in close agree- 
ment with the 170 /im number counts in the small area of 
the "Lockman Hole" (see also Kawara et al. 1998) but their 
90 /im integral counts are significantly above the model. 

The nature and redshift distributions of the ELAIS 
galaxies can test the various models and their hypothesis 
e.g. distinguishing the different galaxy populations on which 
these models are built. This will also help to clarify the origin 
of the differences seen in the number counts of the various 
ISO surveys at different wavelengths. 

The 90 /im luminosity function is presented in Serjeant 
et al. (2004) and Rowan-Robinson et al. (2004) present re- 
sults based on the ELAIS final-band merged catalogue com- 
bining the ISO and ground-based observations in the ELAIS 
fields. 
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